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Abstract 

We analyze the gaps in the excitation spectrum of a Haldane chain 
with single-ion anisotropy in a staggered field. We show that the gap 
along the direction of the field increases at a faster rate than the others, 
while its spectral weight decreases, being transferred to a two-magnon 
continuum. 



In recent times the preparation of a new class of magnetic materials of the 
general formula (R^Yi-a^BaNiOs, where R is one of the magnetic rare-earth 
ions, has arisen a new theoretical interest for the study of integer-spin Haldane 
chains in a staggered magnetic field. The reference compound Y2BaNi05 jl| 
is one of the best examples of a magnetic material with a singlet ground state 
and a spin gap which has been interpreted in the framework of the Haldane 
conjecture Q, according to which integer-spin chains have disordered ground 
state separated by a gap from a degenerate triplet of excitations. In all these 
materials |||, the Ni ions are arranged along linear chains (the a-axis of their 
orthorhombic structure) with a relatively large in-chain antiferromagnetic cou- 
pling (J w 300 K) and negligible inter-chain coupling. When partially or to- 
tally substituted to the yttrium, the magnetic R ions get positioned between 
neighboring Ni chains and are weakly coupled to them. Below a certain Neel 
temperature Tjy (typically: 16 K< Tjy < 80 K) they order antiferromagneti- 
cally. This has the effect of imposing an effective staggered (and commensurate) 
field on the Ni chains, whose intensity can be indirectly controlled by varying 
the temperature below T/v . A phenomenological approach to the study of these 
compounds has been first presented in where the authors propose a simple 
mean field model, in which the bare magnetization curve for the Ni sublattice is 
assumed in the form of the staggered magnetization of a single S=l spin chain, 
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with the effective staggered field proportional to the magnetization of the R 
ions. The so obtained magnetization curve as function of the magnetic field is 
in good agreement with the experimental data. 

The staggered field is expected to partially lift the degeneracy of the Hal- 
dane triplet, leading to different spin gaps in the longitudinal (i.e. parallel to 
the field) and transverse channels. Surprisingly, up to our knowledge there is 
only one experiment || , analyzing spin-polarized inelastic neutron scattering on 
Nd2BaNi05 crystals, which yields a clear evidence for the longitudinal mode. 
In this experiment it is found that at T > T/v = 48 K the value as well as the 
intensity of the transverse and the longitudinal modes are practically equal, as 
in the case of an isolated Haldane chain. To be precise, the resolution of the 
equipment is high enough to resolve two slightly different peaks in the spectral 
density, resulting from a very weak single-ion anisotropy along the chain. Be- 
low T/v, the two excitations clearly separate, with both gaps increasing as the 
temperature is lowered (the effective staggered field is increased) , but with two 
different rates. Not only the longitudinal gap has a much slower increase but 
also its weight seems to be suppressed below T/v, contrary to what happens to 
the transverse mode whose weight is not affected by the magnetic order. 

A theoretical approach to isotropic integer-spin chains in a staggered mag- 
netic field has been first considered in B, starting from the familiar mapping 
[|[ [?| of the Heisenberg model onto a nonlinear cr-model (NLcrM) effective La- 
grangian, that describes the continuum limit, low-energy physics of the Heisen- 
berg chain. What the authors discuss very accurately is actually a related and 
somewhat more phenomenological model in which the strict NLcrM constraint 
is softened, then replacing the original NLcrM with a theory of the Ginzburg- 
Landau type parametrized by an adequate set of adjustable parameters. Some 
recent very accurate numerical studies , based on an extensive Density Matrix 
Renormalization Group (DMRG) analysis of an S=l Heisenberg chain in a stag- 
gered field, have revealed a significant discrepancy with the above mentioned 
theoretical approach, from which it differs in the high field regime and especially 
as far the longitudinal channel concerns. This has led also the authors of ref. 
[D to question the validity of the NLcrM approach. 

This problem has been solved in [H , where it has been shown that an accu- 
rate treatment of the NLcrM does indeed lead to an excellent agreement between 
the analytical and the DMRG results. We will review this treatment in some de- 
tail below, while applying it to the case in which the Hamiltonian contains also 
a single-ion anisotropy term. Let us now recall only its main features. In this 
approach the NLcrM constraint is not softened, but enforced consistently at each 
level of approximation. In the theory there is then just a single free parameter, 
which is assumed to be the (zero-field) Haldane gap and taken from the DMRG 
data. Also, a consistent scheme of calculation of the propagators and their 
analytical structure is developed, which shows remarkable differences between 
the physical properties of the transverse and the longitudinal channels, already 
at the tree-level of a loop expansion. As expected, the transverse propagators 
has only two simple poles so that the theory is purely bosonic, with the spec- 
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tral weight completely exhausted by this magnon excitation. As a consequence, 
the transverse gap At coincides with the Haldane gap. On the contrary the 
longitudinal propagator has a much richer structure. The dynamical structure 
factor has again two well-defined poles corresponding to single-particle excita- 
tions. However they are now weighted by a prefactor which steadily decreases 
as the staggered field is increased. The spectral weight which is lost from the 
poles gets transferred to a multi-particle continuum, which starts from a two- 
magnon threshold u) > 2 At- Thus we are no longer in the range of applicability 
of the so-called single mode approximation (SMA) which establishes a relation 
between the (longitudinal) susceptibility and the gap A^. Once the continuum 
is taken into account, the agreement between the NLcM and the DMRG results 
become remarkable. Let us finally mention that this theoretical scenario has 
also been proposed []l0| to explain recent experimental data on CsNiCL3 
which were not well described by the NLerM on isolated spin chains. 

Keeping in mind the structure of the (RzYi-.r^BaNiOs compounds, our 
microscopic starting point is the Hamiltonian 

H = J2 { J ^ ■ 3+i + K ^ 2 - • (!) 

i 

that describes an antiferromagnetic S=l Heisenberg chain with a single- ion 
anisotropy term along the chain z-axis, in a staggered magnetic field. We as- 
sume K <C J and at the end we will distinguish between the cases in which the 
field is parallel and perpendicular to the anisotropy axis. 

We start analyzing this model by mapping the Hamiltonian ([!]) onto a 
NLctM, by making the well known Haldane ansatz Si = SCti with 



a = (-i)xyi-il + ! (2 ) 

hi representing the slowly- varying local staggered magnetization and ij the local 
generator of angular momentum. In the low-energy continuum limit, the model 
can be described effectively by the action 

S eff = J dx{C a - SH ■ n(x) + R\n(x) ■ z\ 2 } (3) 
C a = ±-{c 2 \d x n\ 2 + \d T n\ 2 ) (4) 



where we have set x = (r, a;), J dx = J Q dx J Q dr and the constants c, <?, R are 

2 
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fixed by the microscopic parameters via the relations c = 2JSa, g — % and 



a. 
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We can now impose the unitary constraint |n(x)| = 1 with the aid of a La- 
grange multiplier field A(x) by writing the partition function Z = TV{exp[— (3H\) 
as the path-integral 
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2tt _ 


exp 


- S eff - i J 



dxA(x)(|n(x)| 2 - 1) 



(5) 
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The isotropic case R = has already been considered in ref. 0, to which 
we refer for all the details of the calculations, that can be easily transposed to 
the case we are interested here. In this letter we will outline only the main steps 
that lead to the calculations of the physical quantities. 

First we need to promote Z to a generating functional Z[J] by replacing 
S e ff with 

S[J] = J dx{C a - S*J(x) • ri(x) + R\n(x) ■ z\ 2 } (6) 

and setting J(x) = H only at the end of the calculations. S[J] being quadratic, 
we can proceed to integrate out the fields n(x) obtaining 

Z[J\<* |[^]exp(-S[A,jj) (7) 



with 



S[\,T\ = Tr{\n{G^)} + 1 -Tr{\n{G^ + R)} (8) 
+ i J dxA(x) - itf 2 J dxdx'J(x) ■ G(x,x')J(x') 

where Tr[A(x, x')] stands for J dxA(x,x). We have denoted with G(x, x') the 
3x3 matrix operator whose inverse G~ 1 (x, x') has components G~^(x, x') (for 
a,@ = 1,2,3) given by 

G~i(x,x') = --^\c 2 d 2 x + d 2 + 2igc\{x) - gcRS a3 ]5( X x') (9) 

Finally we resort to a saddle-point approximation to implement the con- 
straint by fixing A(x) to its mean field value A*(x, J) through the equation 
5S[X, J]/5X(x) = which reads 



1 = 2G n (x,x) + 1 + S 2 [ dydy'J(y) ■ G(y, x)G(x, y') J(y') 

G n (x, x) + R J 

(10) 



In this approximation, we can calculate the connected components of the prop- 
agators as 

P Q/3 fx x') « S 2 S[X*,J] 

r c (x,x)~ sj°(x)8JP{x>) 1 j 

which are therefore given by 

P c ^(x,x') = ^r /3 [G QQ (x,x') + G Qa (x',x)]+*5 2 y'd y d y '[G(x, y )G(y,y') 
+ G(y',y)G(y,x)]^J"(y)|^ (12) 
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The functional derivative 5\(y')/6J@(x) in the right hand side is implicitly 
determined by the integral equation 



J d yj ff(x,y)^^M_) =-2S 2 G QQ (x,x') J dyG aa (x,y)J a (y 



) (13) 



with 



HM = { M\ =4r(X)X0 + ^r 1 1 (x ) x ) r(x ) x') Gri 1 (x ! x) 



5A(x)5A(x')y /iA=A . v ' ; (i? + G n 1 (x,x))2 

+ 45 2 ydydy'J(y)G(y,x')G(x',x)G(x,y')J(y') (14) 

where r(x,x') is the polarization bubble T(x,x') = Gn(x,x')Gn(x,x'). 

At the physical point J = H, one finds that the field A is constant. Setting 
then c 2 £ -2 = —2igcX the following results can be found. 

Going to Fourier transform in q = (q, D, n — 2irn/(3), the sums over the 
frequency can be exactly performed while, in the continuum limit, it is necessary 
to introduce a cut-off A in momentum space. The saddle point condition ( [l0| ) 
then reads 



r A 

1 = 9 - dq 



vV + £- 2 2 
lcoth^/fcyV +g~ 2 + t) 
2 v/g 2 +r 2 +t 



' : coth(i/3<V-7 2 + r 2 ) (15) 



f 4 n 2 a 2 



5 (iJj, + iJ ) — s- + S H 



" z c 2 (t + £,- 2 ) 2 



with i = Rg/c. As in 0], we have chosen to fix the only free parameter of 
our theory, the cut-off A, by requiring the correlation length £~ 2 to coincide 
at H = with the numerical result that can be obtained through a Density 
Matrix Renormalization Group analysis jl2| of the original Hamiltonian (|l|) at 
zero field. 

Then the magnetization can be easily calculated through the formula 

. = 1 « W %' y , (16) 



/3Z<5# c 



Also, assuming translational invariance, eq. ( |l4|) can be solved by going to 
Fourier transform to obtain 

The latter formula can be now inserted in ( |l3| ) in order to calculate the Fourier 
components of the propagators (|l2|). 
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1) Let us first examine the case in which the applied staggered field is parallel 
to the anisotropy axis: H — Hz. In this case m x — m y = 0, while m z as a 
function of H is obtained from (pf). 

One finds that the transverse susceptibility \ T — Xx = Xy is given by 



T _ m _ gS 2 1 



Also, the transverse propagatorP T = P xx = P vv is just a free boson propagator: 

P c T (q) = S 2 G n (q) (19) 

with 

Gll(q)= [n 2 n + c 2 (q 2 +^ 2 )} (20) 

which, when analytically continued to the real axis, has simple poles at u) = 
±ce(q), e(q) = \J q 2 + £~ 2 . Thus the spectral weight function is simply given by 

ImG T c {q, u) = ~ ~ % + *(«))} (21) 

and has therofore the structure required for the applicability of the SMA |l3| . 
In this approximation, \ T — Sgc/(A 2 ^) with the transverse gap coinciding with 
the Haldane gap: 

A T = A = cr 1 (22) 
On the contrary, the longitudinal susceptibility \ L = \z has to be calculated 



as 



L dm T fi , nrr £ ^ 



* = m=* \ 1 + 2H tT^W (23) 



with obtained from eq. ( p^|) 
d£ 2HS 2 g 2 



dH c 2 (t + ^ 2 ) 2 



A 



* v 71 + ( A 2 



(24) 



g Ar 3 AS 2 g 2 H 2 £" 3 1_ ' 



2^(r 2 + i)(r 2 +A 2 + i)^ c2 (^ + r 2 ) 3 

In this case, the longitudinal propagator P L — P zz is instead given by 



5 2 (2 + iy 2 )G 33 (q)r(q) 
(2 + iy 2 )r(q) + 2 ^G 33 (q) 



^(q) = ,:. ' ^t: ; (25) 



G 



with 



G 33 (q) = 

w 2 - 



gc 



[n 2 n + c 2 (q 2 +^ 2 +t)} 
G n 2 (0) 



(26) 



(R + G^(0)f 

Following ||, we can perform an analytic continuation in the variable lo (u 
z = lu + it]) and define: 



gcS 2 



2gc 



's 2 



(27) 



, e(?) = Vc 2 (g 2 +^ 2 +0 



Thus we are led to analyze the structure of the function (omitting the specifi- 
cation of the dependence of the label q) 



G{z) 



Y{z){e 2 - z 2 ) + 5 



(28) 



for < S < 1 and e < 2A . 

Going to the real axis from above ( i] — ► + ), a long but standard calculation 
shows that, for \lu\ < 2Aq, the imaginary part of G{z) is given by 



ImG{u) = j—{5(lo - e L ) - 5(lu + e L )} 
2£l 



(29) 



so that the longitudinal propagators has simple poles on the real axis at the 
points lu = ±€l, with €l fixed by the equation 



ri(ei) 



(30) 



where T{z) = Ti(z) + iY2{z). The prefactor 7 in eq. (|29| ) gives the reduction of 
the quasiparticle weight with respect to the pure bosonic case and is given by 



7 = 



5 dT^Lu 2 ) 



dtu 2 



(31) 



The coefficient 7 is steadily decreasing with increasing H , going to the value 1 
quadratically for H — > 0. 

In the range \oj\ > 2A , one instead finds: 



ImG(u) 



5T 2 (uj) 



(e 2 - lo 2 ) 2 {T 2 {lo)) 2 + (5+ (e 2 - co^T^uj)) 2 



(32) 



As H — > 0, ImG(uj) vanishes and we recover the simple pole structure for both 
the longitudinal and the transverse propagators. The longitudinal pole will also 
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Figure 1: The longitudinal gap (solid line) and the transverse gap At 
(dashed line) as a function of H = Hz for (a) K = -0.02 and (b) K = +0.02. 



survive up to saturation, but with a strongly field-dependent strength: as the 
field increases the spectral weight that is lost from the pole gets transferred to 
the two-magnon continuum (B3) , as it is dictated by the sum rule H : 



J ^ w ImG( w ) = l 



(33) 



As in the transverse case, we can define the longitudinal gap A^ via A^ = 
= 0), thus obtaining 



A| = A 2 T + c 2 t + 



ri(o,A L ) 



with 



Ti(0,A L ) 



g(2 + W 2 



2A T I" v/w 2 - 4A 2 



(34) 



(35) 



Due to the presence in the Hamiltonian (^) of the single-ion anisotropy term, 
the longitudinal and the transverse gaps do not coincide even at H = 0, the 
former being larger (smaller) than the latter for K > (K < 0). The case 
K = —0.02 and K — +0.02 are shown in Fig. 1(a) and Fig. 1(b) respectively, 
where we give A^ = A z and At = A x = A y as function of the magnetic field. 

We notice that for positive K the longitudinal gap is always higher than the 
transverse one. On the contrary, for negative K, the longitudinal gap is lower 
than the transverse one at zero field. As a function of H the former increases 
at a faster rate than the latter, thus resulting in a level crossing at H ~ 0.075, 
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as it is shown in the inset of Fig. 1(a). 



2) Now we examine the case in which the applied staggered field is perpen- 
dicular to the anisotropy axis: H = Hx. We concentrate our attention on the 
structure of the propagators, which now are all different: 

p**(q) = s 2 (2 + ^ 2 )Gii(q)r(g) 

(2 + W 2 )r(q) + ^Gn(q) 
Pf(q) = 5 2 G u (q) (37) 
P-(q) = S 2 G 33 (q) (38) 



where Gu(q), G33(q) and W are given in (|20|) and (J26|) . 

For the transverse channels y, z one immediately sees that the propagators 
are again those of a free bosonic theory to which one can apply the SMA. Thus 
the gaps are simply given by 

A, = Ao^cr 1 (39) 
A z = + cH (40) 

As before, the SMA is no longer applicable in the longitudinal x-direction. 
We find: 

a - = a ' + lMok) <41) 

with 

V (0 A ) g(2 + W2) [ + °° du 1 1 U2) 

For negative K, we find an initial splitting between a lower gap in the z 
direction and a higher one in the x, y directions. As the field is turned on along 
the x-axis, we find a further splitting of A x and Aj,, with the former increasing 
faster. This is shown in Fig. 2(a). On the contrary, for positive K the splitting 
at zero field is between a lower gap in the x, y directions and a higher one in the 
z direction (Fig. 2(b)). As before, the gap A x increases faster than A y and A z 
when the field is turned on along the x-axis. At H ~ 0.06 there is a crossing of 
levels with the gap A^ becoming larger than A z as shown in the inset of Fig. 2(b). 
In both cases, towards saturation the longitudinal gap A x goes to a value about 
twice that of the transverse modes A y , A z , which become approximately equal. 
In addition, though surviving up to saturation, the longitudinal mode has a 
strongly field-dependent strength: as the field increases the spectral weight that 
is lost from the pole gets transferred to the two-magnon continuum. 

In conclusion, we have analyzed the effect of a staggered magnetic field 
on a S — 1 Haldane chain with single-ion anisotropy. We have found that 
the zero-field splitting of the Haldane triplet due to the anisotropy term gets 
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enhanced by the presence of the field. In addition, as the field is increased, 
the longitudinal mode has a spectral weight that decreases quite fast and is 
transferred to a continuum, which has a two-magnon threshold. We remark 
that this seems to be consistent with recent neutron scattering experiments on 
quasi-one-dimensional compounds ||. 
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